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PREFACE 


In  a  previous  report   (AMP  Report  42.11),    the  problem 
of  vertical  impact  of  a  body  with  a  spherical  nose  on  a 
water-surface   was   solved  mathematically  by  Maz  Shlffman  and 
Donald  C.   Spencer  of  the  New  York  University  Group  of   the 
Applied  Mathematics  Panel,      The   solution  was  based  on  the 
assumption   that  the  svirface  of  the  water  is  not   appreciably 
altered  by  the  entry  of  the   body,   an  assumption  which  can 
be   accepted  only  as  a  first  approximation. 

Therefore,   a  further  study  was  undertaken,    taking 
into  account  the  disturbance  of  the  surface  of  the  water. 
The   present  report  presents   the  results  of   these   studies. 
It  is  shown   that  the  disttarbanc©  of  the  surface   leads   to  two 
types  of  corrections   of  the  previous  results.     In   the  first 
place,   the  upward  rise   of   the   svirface  of  the  water  wets  a 
larger  portion  of   the  entering  body,   thereby  increasing   the 
resistance  of  the  water  to  the  entering  body.     Secondly,    the 
rising  surface  relaxes  somewhat   the  restraint  Imposed  on  the 
body,   causing  a  decrease  in  the   Impact  force.      These   two 
corrections,  which  partly  counterbalance  one  another,   are 
analyzed  here  by  using  a  combination  of  theory  and  information 
derived  from  experiment.      The  final  result,  which  deviates 
somewhat  from  the  first  approximation,   is   compared  with 
experiment  and  shown  to  be  in   satisfactory  agreement  with  it. 

With  the  experience  gained  and   the   theory  developed  for 
vertical  entry,   it  is  hoped   that  the   Impact  forces  during 
oblique   entry  will  be    amenable    to   treatment.      This   question 
is  now  being  considered  by  the   same  authors. 

Attention  should   be  drawn  to  the  work  of  the   Harvard 
Group  of   the  Applied  Mathematics  Panel  concerning  air-water 
entry  as  a  whole,   and  to  the  activities  of   the  Underwater 
Ballistics  Committee. 

The  writers  were  greatly  helpied  by  other  members  of  the 
New  York  University  Group,   in  particular  by  Mr.  M.    H.   Shamos, 
who  made  an  important  contribution  by  excellent  flash  photo- 
graphs he   took  of   spheres  entering  water.      These   photographs 
fxirnished  much  qualitative  and  numerical  information.      Thanks 
are   also  duo    to  Dr.    Patrick  Hurley  of   the  Morris  Dam  Group  of 
the  California  Institute  of  Technology  for  permission   to  use 
the   photographs   exhibited   in  Plate    5;      to  Dr.   R.  M.   Davles  of 
the  Engineering  Laboratory,   Cambridge,   England   for  Plate   6; 
and  to  Dr.  E.   Newton  Harvey  of  Princeton  University  for  Plate  ' 

R.   Courant 
Director,   Contract  OEMsr-945 
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THE  FORCE  OF  IMPACT  ON  A  SPHERE  STRIKING  A  WATER  SURFACE 
Second  Approximation 

Section  1.  Introduction  and  Summary 

The  entry  of  a  solid  body  from  air  Into  water  may 
take  place  In  two  ways:   (1)  If  the  speed  Is  small  and  If 
the  siarface  of  the  body  Is  smooth,  the  entry  occiars  with- 
out the  formation  of  a  cavity  (smooth  entry);   (11)  If  the 
speed  Is  great  or  If  the  surface  of  the  body  Is  rough  or 
has  an  Irregular  shape,  the  flow  detaches  from  the  body 
when  It  has  penetrated  a  short  distance  below  the  Initial 
surface  level  (rough  entry).  In  either  case  there  Is  an 
Initial  period  following  contact  with  the  surface  during 
which  the  body  experiences  the  greatest  deceleration  and 
which  Is  of  such  short  duration  that  cavitation  has  not 
yet  developed*  The  determination  of  the  forces  acting 
during  this  Impact  stage  of  the  entry  Is  Important  In  the 
case  of  projectiles  shot  or  dropped  Into  water,  not  only 
because  of  possible  damage  to  the  mechanism  and  nose  but 
also  because  the  In^julslve  forces  and  torques  created  by 
the  Impact  Influence  the  underwater  trajectory  of  the  pro- 
jectile. 

We  suppose  that  the  nose  of  the  projectile  Is  spheri- 
cal In  shape  and  that  the  Initial  Inipact  with  the  surface  Is 
vertical.  We  may  then  visualize  the  mass  of  the  projectile 
as  concentrated  Into  a  sphere  of  the  same  radius  as  the  nose. 
In  a  previous  report  [15]  we  obtained  a  first  approximation 
to  the  force  of  Impact  on  a  sphere  entering  vertically  by 
disregarding  the  rise  of  the  surface;  in  the  present  report 
we  obtain  a  second  approximation  by  estimating  the  effects 
of  the  surface  motion  and  we  compare  our  theoretical  estimates 
with  experiment. 


*  Numbers  In  brackets  [  ]  refer  to  the  bibliography  at  the  end 
of  this  report. 
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The  resultant  upward  force  P  acting  on   the  sphere 
In  the  Impact  stage  can  be  expressed  In  terras   of  a  dlmen- 
slonless   Impact-drag  coefficient  C  .      Let  A  be   the  radius 
of  the    sphere,  U     its   Initial  velocity  upon  Inpact,  and 
f>  the   density  of  water.      Then  the   Impact   force  may  be 
conveniently  expressed  as 

where  C     is    the   "irapact-drag  coefficient."     The   coefficient 
C     depends  on  the    time;     it  is   zero  when  the   sphere  first 
strikes   the  surface,  rapidly  Increases   to  a  maximum  value, 
and  then  decreases. 


The   Impact-drag  coefficient  C     may  be   expressed  In 
terms   of   the  "virtual  mass"   of   the  fluid.      If  U  Is   the 
velocity  of  the  sphere  at  time  t,  we  define  the  virtual 
mass   of   the  fluid  to  be  a   quantity  M  such  that  MU  Is  the 
total  vertical  Impulse  communicated  to  the  fluid  during 
the  time  Interval  t;      In  other  words  M  Is  defined  by  the 
eqxiatlon 

P  =^   (MU)    . 

Let  B  denote  the  depth  of  penetration  of  the  bottom  of  the 
sphere  below  the  Initial  surface  level,  and  let  b  =  B/A. 

CONFIDENTIAL 


CONFIDENTIAL 


Introducing  the  additional  dlmenslonless   quantity  m  =  M/^trp  A  , 
we  obtain  In  Section  3  the  formula 

dm 

M 


°'  °  (1  .  §)' 

where    (T   is  the   "specific  gravity"   of   the  incoming  projectile 
considered  as    though  the  whole  mass  were  concentrated  into 
a  sphere  having  the  same  radius  as   the   nose, 

J-    _  mass  of  projectile 
-ir^A 

Our  final  estimate   of   C     as  a   function  of  b   is  plotted 

in  Graph  1  for  various  values   of   c  •    (All  graphs  are   collected 

together  at  the  end  of   the  report.     For   the   sake  of  convenience. 

Graph  1  is  also  inserted  on  the   next  page.)     For  values   of  <r 

ranging  from  1  to  oo  ,    the  impact  drag  coefficient   C     rises  very 

rapidly  to  a  peak  of  about   1,   the  peak  occ\irring  when  the  sphere 

has  penetrated  to  a   distance  of  from  0*1  to  0.2  of  a  radius 

below   the   initial   surface   level.      The   coefficient   C     then  de- 

P 
cllnes  more  slowly,   until  at  about  0.7   of  a  radius   of  entry 

its  value   is   equal   to  the   stationary  cavity  drag  value   of  0.3 

(for   description  of  cavity  drag,    see   [2]   aiad  [21]). 

For  smaller  values   of   (T  ,   this   general  behavior  of  C 
is   retained  but   is   shifted  downward  and  to  the   left. 

The  question  arises  whether,   in  rough  entry,    the   cavity 
formation  which  takes  place  will  interfere  with  the   above  re- 
sults and  partly  vitiate   them.      It  is  known  from  both  experimen- 
tal and  theoretical  grounds   that  the  cavity  begins   to  form  only 
after  the  sphere  has   submerged  a  dls  tance  somewhere  between  a 
quarter   of  a   radius  and  a   radius.      Thus    the  behavior   of  C     de- 
picted in  Graph  1  still  applies    the   rapid  rise   of  C     to  a 
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maximum  of  about  1,  and  the  slower  decline  to  the  cavity 
drag  value  of  0.3»   The  brief  ln5)act  stage  In  rougji  entry 
Is  the  same  as  In  smooth  entry. 

The  dlmenslonless  virtual  mass  m  defined  by  m  =  }A/—vpA 
Is  plotted  as  a  function  of  b  In  Graph  3»  For  purposes  of  com- 
parison we  have  plotted  In  the  same  graph  our  first  approxima- 
tion for  m,  taken  from  the  previous  report  [15].   Likewise, 
In  Graph  2,  our  corrected  estimate  of  C  In  the  case  cr  a=  oo 
Is  compared  with  our  first  approximation  obtained  in  [15].  In 
the  corrected  estimate,  the  maximum  of  Cp  is  somewhat  larger 
and  is  displaced  to  the  left. 

The  first  approximation  of  the  impact-drag  coefficient 
was  derived  In  our  report  [ 15]  under  the  assumption  that  the 
surface  does  not  rise,  the  surface  remaining  plane  and  at  its 
initial  level.  Under  these  circumstances  the  flow  at  each 
instant  t  depends  only  on  the  depth  of  penetration  below  the 
initial  level  and  on  the  velocity  of  the  sphere  at  that  Instant 
t  and  is  otherwise  Independent  of  the  past  history  of  the  entry; 
for  this  reason  the  flow  might  be  called  "quasi-steady." 

In  the  present  report  a  second  approximation  is  obtained 
which  takes  into  account  the  past  history  of  the  phenomenon. 
The  corrections  are  of  two  sorts.  On  the  one  hand  the  water 
rises  and  wets  a  larger  portion  of  the  spherical  surface  than 
it  would  If  the  surface  did  not  rise.  The  fluid  therefore 
exerts  pressure  over  a  larger  area  and  Increases  the  total 
impact  force.  This  will  be  called  the  "wetting  correction." 
On  the  other  hand,  the  free  surface  rises  and  relaxes  to  a 
certain  extent  the  restraint  Imposed  on  the  sphere.  This  will 
be  called  the  "surface  correction."  These  two  corrections 
are  discussed  in  Sections  6  and  7  respectively,  and  both  are 
combined  in  Section  8  to  give  the  final  estimate  of  the  impact- 
drag  coefficient  Cp. 

In  Section  9  this  final  estimate  of  CL  Is  compared  with 
experimental  measurements  of  the  impact  force  on  a  sphere. 
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Coii5)arl3on  la  made  with  measurements  taken  by  S.  Watanabe'a 
[17]   In  1930,  who  used  a  piezoelectric  gauge  to  determine 
the  force;     with  rough  experiments  of  R»  W»  Blundell  [3]   made 
In  1937;     and  with  recent  measurements  of  E.  G-   Richardson  [13]  • 

Watanabe's  experiments  were  very  carefully  performed 
with  acovirate  Instruments   and  verified  the  fact   that  the   Impact 
force  varies  as  the  square  of  the   Initial  velocity.     The  value 
of  the  specific  gravity    <r   of  the    sphere  In  Watanabe's  experi- 
ment was  as   low  as  ©•116»      The  theoretical  C_   curve  correspond- 
ing to  this  value   of    (T    Is   compared  with  Watanabe's   experimental 
curve  In  Graph  10.      Note  that  the   range  of  b  Is   small,  from  0 
to  •03.      The  agreement  Is  remarkably  good  for  very  small  values 
of  b,  when  C     Is  Increasing  rapidly.      In  fact,   near  b  =  0     the 


theory  gives 


C     =  6.62  b^/^   -  agb  +  •, 


where  the  constant  0.2  ^^   ^^^  accurately  determined  by  the  theory. 
The  first  term  6.62b  V^    of   this   expansion  Is   verified  completely 
by  Watanabe's  data.      The  deviation  In  the  remaining  portion  of 
the  curve  Is  dae   to  the  second  coefficient  ag. 

In  Blundell 's  experiment,    <r   ranged  from  .2  and  .3  approxi- 
mately.     The   experimental  curves,    transcribed  to   the  dlmenslon- 
less  form  of  a  drag  coefficient,  are  shown  In  Graph  11(a),    and 
the  corresponding  theoretical  curves  are  drawn  In  Graph  11(b). 
An  Idea  of  the  rather  large  experimental  error  may  be  seen  by 
comparing  curves   (B)   and  (C)  which  should  be  identical,   and  from 
the  fact  that  curve   (D)   should  be   tl»   highest   (since    <r   Is   the 
largest).     The  experiments  at  least  give  an  idea  of  the   order  of 
magnitude  of   the  In^jact  fijrces  and  of  their  duration,  and  the 
agreement  with  theory  Is  satisfactory.      In  fact,    curve    (D)    even 
agrees  numerically  with  theory. 

In  Richardson's  experiment,   cr  =  .16.      The   corresponding 
theoretical  curve  is   compared  with  Richardson's   in  Graph  12, 
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and  the   agreement   Is   fair.      Again,   Richardson's   experiment 
should  be  considered  as   yielding  only  the   order   of  magnitude 
of   Cp. 

In  each   of    these  experiments,    the  value   of    C    Is   small* 
It  would  be  desirable   to  carry  out  experiments  with  larger 
values   of    (T  ,   using  accurate   Instruments    (such  as  piezoelectric 
gauges). 

In  addition  to  the   impact   force,    it   is   of  interest   to 
know  the   average  pressure  p   acting  over   the  portion  of   the 
sphere   in  contact  with  the  water.      This  pressure   can  be  written 
In  the  form 

P  =1/°^   •    I>p 

where  D  is  a  dimenslonless  coefficient.  The  dependence  of 
Dp  on  b  =  B/A  Is  shown  in  Graph  4  (for  the  case  (t  =  ©o  ).   It 
Is  noteworthy  that  the  value  of  D_  at  the  initial  contact  of 
the  sphere  with  the  water  is  infinite.   Actually,  because  of 
the  oompressibility  of  the  water  and  because  of  the  fact  that 
the  compression  wave  is  Initially  a  plane  wave,  the  pressure  p 
during  this  phase  Is  approximately  cUq  where  c  is  the  speed 
of  sound  In  water.  This  would  give  a  value  of  Dr,  approximately 
equal  to  -^^  •  The  Intersection  of  the  straight  line  Dp  =  -^ 
with  Graph  4  gives  an  indication  of  the  time  Interval  in  which 
con5)res3lblllty  plays  a  role.   Thus,  the  effect  of  compressibility 

♦  Using  the  asymptotic  formula 

for  small  b,  and  equating  this  to  •—  ,  one  obtains 

b  = 

as  an  estimate  of  the  relative  depth  of  penetration  in  which 
con5)res3lblllty  plays  a  role.  Even  if  the  initial  entry 
velocity  Uq  Is  as  high  as  500  ft/sec  this  yields  b  =  .01,  an 
order  of  magnitude  which  Is  negligible. 

This  Idea  was  also  suggested  by  M.  Glmprich. 
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l3  to  reduce  the  value  of  Du  and  therefore  of  C  at  the 
very  beginning. 

Concerning  the  actual  distribution  of  pressure  over 
the  wetted  portion  of  the  sphere,  we  remark  that  the  pressure 
is  not  greatest  at  the  lowest  point  of  the  sphere  (stagnation 
point).  This  fact  arises  from  the  non -stationary  character 
of  the  flow,  and  the  usual  conceptions  based  on  stationary 
flow  do  not  apply.  On  the  contrary,  the  pressure  is  a  mini- 
mum at  the  lowest  point  of  the  sphere  and  has  a  very  sharp 
and  high  maximum  near  the  edge  of  the  wetted  portion  of  the 
sphere*  This  fact  was  derived  theoretically  by  Wagner  [16] 
in  the  two-dimensional  case  of  the  entry  of  a  very  flat  wedge- 
However,  an  estimate  of  the  depth  at  which  stationary  flow 
approximately  sets  in  can  be  obtained  by  determining  when  C 
Is  about  0.5.  This  occurs  for  b  =  .7  approximately. 

Section  2.  Photographs  of  Entry 

The  interesting  sequence  of  events  produced  by  entry 
of  a  sphere  from  air  into  water  was  first  photographed  and 
described  by  Worthington  [19].  Plate  1  is  a  reproduction 
of  some  of  his  photographs  and  shows  several  stages  in  the 
smooth  entry  of  a  sphere*  A  more  complete  sequence,  which 
shows  the  splash  in  greater  detail,  is  given  in  Plates  2,3. 
As  may  be  observed  from  these  two  sets  of  pictures  the  phenom- 
ena accompanying  smooth  entry  are  closely  duplicated  in  differ- 
ent experiments.  Even  the  Jets  produced  as  the  water  envelopes 
the  spheres  differ  little  in  appearance. 

In  order  to  show  more  clearly  the  splash  during  in5)act, 
an  enlarged  photograph  is  reproduced  in  Plate  4.  We  observe 
that  the  splash  is  largely  concentrated  in  the  immediate 
neighborhood  of  the  sphere  and  is  composed  of  a  relatively 
thick  base  from  which  a  thin  sheath  rises  up  along  the  sides 
of  the  sphere.  A  few  drops  which  have  detached  themselves  may 
be  seen  at  the  top  of  the  sheath. 
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Because  of  the  speed  with  which  projectiles  are 
usually  launched  into  water  their  entry  Is  rough.   In  fact, 
the  entry  of  a  bomb,  mine,  or  torpedo  takes  place  In  three 
stages:   (l)  the  Initial  Impact  stage  which  begins  when  the 
projectile  touches  the  surface  and  ends  when  cavitation 
developes;   (11)  motion  through  the  water  with  a  trailing 
cavity  extending  to  the  siirface;   (111)  motion  after  the 
cavity  Is  closed.   The  closure  of  the  cavity  takes  place 
either  at  the  svirface  or  below  the  surface,  but  In  either 
case  the  projectile  carries  with  It  a  pulsating  air  bubble 
which  gradually  disappears. 

The  Initial  stage  of  the  rough  entry  of  a  sphere  Is 
shown  In  Plate  5.   The  sheath,  which  persists  alongside  the 
sphere  In  the  case  of  smooth  entry,  detaches  In  rough  entry 
and  may  break  up  Into  spray.   However,  the  thick  base  of  the 
splash  remains  In  contact  with  the  sphere  and  Is  the  same  as 
in  smooth  entry.   So  far  as  estimation  of  Impact  force  Is 
concerned,  the  effect  of  the  sheath  Is  negligible  compared  to 
that  of  the  thick  base. 

The  further  stages  of  rough  entry  are  shown  in  Plate  6. 
Unfortunately  only  the  first  photograph  shows  the  sphere  in  the 
impact  phase  and  here  the  splash  cannot  be  seen  very  clearly. 
However,  the  sheath  can  be  seen  as  a  light  band  surroxinding  the 
sphere  above  the  surface. 

For  the  usual  projectile  entering  water  the  speed  does 
not  exceed  five  or  ten  percent  of  sound  speed  in  water.   It  is 
therefore  reasonable  to  neglect  compressibility  effects.   For 
the  sake  of  general  interest,  we  include  Plate  7,  which  shows 
a  sphere  entering  at  roughly  half  the  speed  of  sound  in  water. 
Here  the  effect  of  the  compressibility  is  of  course  inportajit. 

Section  5.  Virtual  Mass 

The  impact  drag  coefficient  C  will  be  obtained  in  terms 
of  the  "virtual  mass"  M  of  the  fluid.   Let  U  be  the  downward 
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velocity  of  the  sphere  at  any  Instant  t,  and  let  P  be  the 
resultant  upward  Impact  force  exerted  on  the  sphere.  Then 
the  sphere  exerts  a  downward  force  P  on  the  fluid  and  Im- 
parts momentum  to  the  fluid.  The  "virtual  mass"  M  of  the 
fluid  at  a  given  Instant  t  Is  defined  by 

(3«1)  P  a  ^(MU)   or   MU  «y  P  dt   • 


Let  M^  be  the  mass  of  the  Incoming  projectile.  Neg- 
lecting gravity  compared  to  the  Impact  force  P,  we  have 


(3.2)  .P=^(M^TJ) 


Adding  (3.1)   and  (3.2)  we  obtain  ^  [(m  +  M^)  U]    =  o  and  so 
Integrating, 


(3.3)  (M  +  Mq)U  =  M^Ujj       or       U 


1  *  M/M, 


whea?e  TJ  Is  the  velocity  of  the  Incoming  sphere  at  the  Initial 
instant  of  contact  with  the  fluid. 

Let  B  =  B(t)  be  the  distance  at  time  t  of  the  lowest 
point  of  the  sphere  from  the  Initial  surface  level;  the  velocity 
U  of  the  sphere  Is  then  B  =  -^  •  By  differentiating  (3.3)  with 

*  Thus  MU  is  the  vertical  Impulse  contributed  to  the  fluid  by 
the  moving  sphere.   It  Is  also  possible  to  introduce  a  virtual 
mass  based  on  the  kinetic  energy  of  the  fluid,  but  for  the  case 
of  iinpact  it  differs  from  the  virtual  mass  defined  above.  This 
Is  to  be  expected  from  general  considerations  of  inelastic  im- 
pact. We  shall  use  throughout  this  report  the  definition 
(3.1)  of  the  virtual  mass,  based  on  momentum;   Justification 
of  this  choice  is  given  In  the  Appendix. 
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respect  to  t  and  using  (3«2),  It  follows  that 
(3,4) 


P  =  dM/dB  ^2  ^   dM/dB   ^   ^ 


Let  A  denote  the  radius  of  the  incoming  sphere,  and  let  p   be 
the  density  of  water.   Introducing  the  dimensionless  quantitlt 


(3.5) 


i^ 


,  b 


B 


we  may  write 
(3.6) 


P  =  i  pT^      .  Tf  A^  .  C_ 


2  '  o 


where 
(3.7) 


dm 
n  = Si— 

P         M   3 

(1  -i-r:  ) 
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For  incoming  projectiles  of  large  mass  compared  to 
the  mass  of  the  water  displaced  by  the  hemispherical  nose, 
^  may  be  neglected  coinpared  to  1.  Then  the  impact  drag 
cSefficlent  C  is  merely  '^.      Also,  from  (3»3),  the  change 
in  velocity  U  in  the  impact  stage  is  small. 

However,  if  M^  is  not  large  compared  to  M,  the  denom- 
inator in  (3«2)  must  be  retained.  This  is  the  case  in  the 
experiments  which  will  be  cited  in  Section  9,  and  it  is  con- 
venient to  put  the  denominator  in  a  slightly  different  form. 
Let  (T  be  the  specific  gravity  of  the  incoming  projectile 
considered  as  though  the  whole  mass  were  concentrated  into  a 
sphere  (if  the  incoming  projectile  is  not  already  a  sphere), 
i.e., 

(3.8)  -     "' 


I-m' 


Using  (3.5),  equation  (3.7)  maybe  written 

Am 
(3.9)  C   =  ^ . 


P 


<i*f.»'' 


The  determination  of  the   impact-drag  coefficient  C     requires 
merely  the  determination  of  the   dlmensionless  virtual  mass  m. 

Section  4.  Hydrodynamlc  Formulation  of  the  Problem. 

We  shall  neglect  the  effects  of  compressibility  and  of 
the  viscosity  of  the  fluid.   Since  the  fluid  is  at  rest  before 
the  sphere  strikes  it,  the  flow  of  the  fluid  at  each  instant 
of  time  t  is  irrotational  and  can  be  described  by  a  velocity 
potential  1/ •  This  is  a  function  v(x,  y,  z,  t)  defined  over 
the  region  of  space  occupied  by  the  fluid  such  that 

-  grad  1/ 
is  the  velocity  of  the  fluid  at  the  Instant  t  and  the  position 
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X,  7,  z.   The  velocity  potential,  on  account  of  the  Incompress' 
iblllty  of  the  fluid,  satisfies  the  potential  equation 


■V^^c  *  ^yy  *  ^zz  =  0  ' 


where  subscripts  mean  differentiation.   The  potential  Is 
normalized  by  setting  V  =  0  at  oo  ).   The  pressure  In  the  fluid 
Is  connected  with  the  potential  function  by  the  Bernoulli 
equation 


(4.1) 


£*i(grad^)^.||r=0 


where  p  is  the  pressure  In  excess  of  hydrostatic  pressure* 


'/J  -/o/a^e 
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To  obtain  the  Impact  force  exerted  on  the  sphere, 
we  shall  use  the  principle  of  conservation  of  momentum. 
The  total  momentum  of  the  fluid  Is  directed  vertically- 
downward  because  of  the  rotational  symmetry  and  Is  equal 
to 

(4.2)    //f/°  If  <3xdyd2  =^ffy/dx&j   -/of/'v^dxdy  • 

Here  W  Is  the  wetted  portion  of  the  sphere,  S  Is  the  free 
surface,  and  the  last  term  In  (4.2)  means  the  limit  of  the 
Integral yy'^dxdy  extended  over  the  submerged  portion  H  of  a 
large  spherical  surface  the  radius  of  which  tends  to  Infinity. 

The  external  forces  acting  on  the  fluid  contained  in 
the  membrane  H  are;   (1)  the  downward  force  of  magnitude  P 
exerted  by  the  sphere  and  (11)  the  upward  force  exerted  on  H 
by  the  pressure  of  the  fluid  external  to  it.   (The  resultant 
force  due  to  the  constant  atmospheric  pressure  is  zero.  The 
only  effect  of  gravity  is  the  buoyancy  which  is  negligible 
compared  to  the  Inipact  force  P.)  The  force  (11)  is  equal  to 
the  time  rate  of  change  of  the  last  term  in  (4,2).  For,  by 
Bernoulli's  law  (4#1),  the  force  (11)  differs  from  the  time 
rate  of  change  of  the  last  term  by  a  term  of  the  form 


// 


igcaA-^)^   dxdy 


// 


which  approaches  zero  as  the  radius  of  H  becomes  infinite. 

1  2 

(For,  at  oo  ,  1^  behaves  at  worst  like  -i  , therefore  igr&dyr) 

T  1 

like  -=4  and  the  whole  integral  like— 3)  •  Likewise  the 

moment\im  flux  through  the  surface  H  consists  of  integrals  of 
squares  of  velocities,  and  these  approach  zero  as  the  radius 
of  H  becomes  infinite.  Since  the  resultant  of  all  the  ex- 
ternal forces  is  equal  to  the  time  rate  of  change  of  the 
momentum  (4.2),  we  therefore  have 

(4.3)  P  =  ^  (/>//v^dxdy)  . 
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The  formula  (4»3)  can  also  be  Inferred  directly  from 
the  knovra  physical  Interpretation  of  the  potential  function 
11/   ,  namely  that  yo^  Is  the  Impulsive  pressure  required  to 
produce  the  flow. 

It  Is  convenient  to  exhibit  separately  the  effect  of 
the  velocity  B  of  the  sphere  by  setting 

(4.4)  -^(x,  y,  z,  t)  =  B(l)(x,  y,  z,  t)  , 

so  that  ^  is  the  Instantaneous  potential  function  correspond- 
ing to  the  parameter  B  and  to  unity  as  velocity  of  the  sphere. 
Then  (4.3)  can  be  written 

P-^  (MB) 

where 

(4.5)  M  =  ^fJc^dxAj     . 

W+S 
Mils  formula  expresses  the  virtual  mass  of  the  fluid  In  terms 
of  the  potential  function  <^  • 

Section  5.  The  First  Approximation  for  a  Sphere 

Consider  the  flow  produced  in  the  fluid  as  a  result 
of  the  entry  of  the  sphere.  As  the  sphere  penetrates  the 
surface,  the  free  surface  no  longer  remains  plane  but  rises 
slightly  and  forms  spiirts  running  up  the  sides  of  the  sphere. 
Although  the  exact  shape  of  the  surface  is  difficult  to  deter- 
mine theoretically,  some  estimates  will  be  presented  in  Sec- 
tion 6.  For  the  purpose  of  computing  the  virtual  mass  and  the 
impact  force,  we  shall  make  a  plausible  first  approximation 
in  the  present  section  by  disregarding  the  rise  of  the  surface. 
More  exactly,  we  shall  neglect  the  rise  of  the  sxirfaoe  and 
also  the  squares  of  velocities  on  the  stirface* 
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To  determine   the  potential  function  T/^  (or  ^  ),    It 
is   necessary  to  know  boundary  condition  for    -y/^  •   The   exact 
boundary   condition  on  the  free  surface  S  of  the  fluid  is 


(presstire  is  atmospheric),  which  by  Bernoulli's  equation 
(4»1)    can  be  expressed  in  terras   of  the  potential  function 
F  ^7 


(5.1)  1^:=  1   (gradl^)' 


Since  this  exact  boundary  condition  is  non-linear  and  diffi- 
cult to  apply,  we  neglect  the  rise  of  the  surface  S  and 
squares  of  velocities  on  S  as  a  first  approximation.   Thus, 
neglecting  (gradl^)  in  (5»1)  we  have 


§1=  ° 


or, 

original  level. 


(5»2)  y/  =  0  on  the  initial  plane  z  =  0  (since  1^=  0  at  oo  )• 

The  boundary  condition  (5«2)  states  that  the  inpulsive 
pressure  /o-^  is  zero  on  the  surface.  This  condition  would  be 
correct  if  the  impact  took  place  completely  impulsively* 
Actually  the  time  interval  is  very  short,  and  so  (5.2)  is  a 
good  first  approximation. 

It  is  more  convenient  to  use  the  potential  function 
^      for  unit  velocity,  defined  by  (4«4).   The  potential 
function  (|)  (depending  on  B)  satisfies  the  following  conditions 
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(a)   On  the  sphere,    the  negative   of   the   normal  derivative 
of    y  =  B^)     is  equal  to  the  coinponent  of  the  velocity  of 
the   sphere  along  the  nornial,   or 


-  ii  = 


an 


(b)  On  the  initial  water  level, 

(()  =  0  . 

There  is  a  unique  potential  function  satisfying  these  con- 
ditions. 

The  flow  arising  from  the  above  boundary  conditions 
has  been  constructed  explicitly  in  the  preceding  report  [15]. 
Because  of  condition  (b),  the  potential  flow  can  be  extended 
to  the  infinite  space  by  reflection  on  the  surface  and  the 
flow  at  each  instant  corresponds  to  the  steady  flow  about  a 
symmetrical  lens  formed  by  the  intersection  of  two  equal 


CONFIDENTIAL 


CONFIDENTIAL  18 


spheres.   A  photograpli  from  Worthlngton  [  l9]  Is  shown  below. 
This  photograph  Is  a  view  of  the  entering  sphere  from  slightly 


"below  the  surface.  The  stirface  reflection  creates  the  upper 
side  of  the  lens.   (Of  co\irse,  the  size  and  shape  of  the  lens 
varies  with  time.) 

The  preceding  report  [15]  supplies  not  only  the  potential 
function  (j)  but  also  the  corresponding  virtual  mass,  which  we 
shall  henceforth  denote  by  M^^,  the  subscript  L  referring  to  a 
lens.  The  results  of  [15]  which  we  use  here  may  be  summarized 
as  follows. 


(i)  The  diraensionless  virtual  mass  m--   «  ■  **         is  a  function 


of  b  alone  (where  b  =  B/A),  and  will  be  written  as  Mj^(b).   Its 

dependence  on  b  is  drawn  in  Graph  3.  The  impact-drag  coefficient 

dra-r  (b ) 

— jg —  depends  on  b  in  the  manner  indicated  in  Graph  2.   The 

maximum  of  ^L  is  .95  and  it  occurs  for  b  =  •24. 
do 
When  b  is  small,  the  spherical  segment  beneath  the  initial 

level  of  the  fluid  is  nearly  a  disc.  By  a  known  formula  for  the 

virtual  mass  M  of  a  disc  (see  [8],  p.  130)  we  therefore  have 


!/"»' 
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for  B  small,  where  c  Is  the  radius  of  the  circle  of  Inter- 
section of  the  sphere  with  the  original  level. 
This  gives  o   -K/o 


(5.3)        m^^b)  =  1^  2^2  1,3/2  ,  g^^Q  ^3/2 
for  small  b«   The  more  exact  lens  flow  gives 

'   mj^(b)  =  2.40  b^/^  -  1.15  b^  +  ., 


(5.4) 


^   =  3.60  b-"-/^  -  2.30  b 
db 


(11)  The  variation  of  the  virtual  mass  with  the  depth  of 
penetration  can  be  better  visualized  by  setting 

(5.5)  M^  =  kj^yo  c'^  , 

where  the  dlmenslonless  quantity  k,  depends  on  b.   The  coefficient 
kj.  Is  drawn  In  Graph  5.  When  b  is  small  the  spherical  segment 
Is  nearly  a  disc  so  that 

4 

ICj^-o-   as   b-->0. 

When  b  =  1,  the  spherical  segment  Is  a  hemisphere  and  It  Is  then 
well  known  that  k^  =  ^  (see  [8],  p.  116).  The  quantity  k^  de- 
creases rapidly  from  the  value  —  at  b  =  0  and  then  levels  off, 
even  swinging  upward  slightly  near  b  =  1.  This  behavior  of  k^ 
Is  plausible  on  Intuitive  grounds.  (The  final  coefficient  k, 
which  Includes  the  various  corrections  made  in  this  report,  is 
also  drawn  in  Graph  5. ) 
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It  Is  interesting  to  note  that  for  the  case  of  a  cone 
Instead  of  a  sphere  the  quantity  k  is  a  constant  independent 
of  the  depth  of  penetration  of  the  cone.   This  results  from 
the  fact  that  for  a  cone  the  flow  patterns  at  different  times 
are  geometrically  similar. 

To  obtain  a  second  approximation  to  the  virtual  mass  and 
to  the  impact  force,  it  is  nex3e3sary  to  estimate  the  effects  of 
the  motion  of  the  surface*  On  intuitive  grounds,  corrections 
due  to  the  motion  of  the  surface  can  be  divided  into  two  parts. 
On  the  one  hand  the  water  rises  and  wets  a  larger  portion  of 
the  sphere,  the  sphere  now  exerting  pressure  on  the  fluid  over 
a  larger  area  and  causing  an  increase  in  the  total  impact  force. 
This  will  be  called  the  "wetting"  correction.  On  the  other  hand, 
the  free  surface  of  the  fluid  rises  and  relaxes  somewhat  the 
restraint  imposed  on  the  sphere.   This  will  be  called  the  "sur- 
face" correction. 

These  two  corrections  appear  in  formula  (4, 5)  when  written 
as 

(5,5)         M  =  /oJT^6xd.j   +p/J  (|)dxdy  • 

IV  S 

The  second  Integral  is  not  zero  since  the  condition  (5.1)  re- 
places if)  =  0.   (In  fact,  ^   will  turn  out  to  be  negative  on  S.  ) 

Section  6.   The  Wetting  Correction 

The  greatest  disturbance  of  the  surface  of  the  fluid  occurs 
alongside  the  sphere.  Here  a  thin  sheath  of  the  fluid  rises  more 
rapidly  than  the  rest,  and  may  detach  from  the  sphere  (see  Plate  4). 
This  sheath  because  of  its  thinness  contributes  very  little  to 
the  impact  force  (the  pressure  of  the  fluid  in  the  sheath  is  practi- 
cally atmospheric).  The  sphere  exerts  pressure  on  the  fluid  in 
the  main  over  the  portion  of  the  sphere  obtained  by  cutting  off 
the  thin  sheath  from  the  thick  "base"  of  the  risen  fluid.   This 
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portion  of  the  sphere  is  Indicated  by  the  points  QQ  in  the 
diagram,  and  the  separation  of  thin  sheath  from  the  thick 
base  Is  Indicated  by  the  dotted  lines  QQ'.  The  flow  of  the 
fluid  Is  due  to  the  portion  QQ  of  the  sphere,  and  Is  approxi- 
mately the  flow  about  a  lens  consisting  of  the  spherical  seg- 
ment QQ  and  a  horizontal  plane  free  surface  QQ'oo  •   (This 
may  In  fact  be  taken  as  the  mathematical  definition  of  the 
point  of  demarcation  Q  between  sheath  and  base. )   In  other 
words,  the  "base"  below  QQ»  Is  so  thick  that  Its  Influence  on 
the  sphere  Is  practically  the  same  as  If  It  extended  to  oo  • 
It  will  be  shown  later  In  this  section  that  the  above  discussion 
is  exactly  verified,  and  the  definition  of  Q  Justified,  in  the 
asymptotic  case  as  B  -^  0  • 

*  For  the  method  used  in  the  practical  determination  of  this 
point  of  demarcation  Q  from  photographic  data,  see  Plate  4 
and  the  accompanying  discussion  on  page  33  • 
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Henceforth  It  will  be  the  spherical  segment  QQ  which 
will  "be  called  the  "wetted"  portion  W,  or  the  "effective" 
portion,  of  the  sphere.  Let  B,  be  the  depth  of  this  effective 
portion,  and  c,  the  radius  of  Its  circular  base  (see  diagram). 
The  quantity  B,  will  be  called  the  "effective  depth"  of  penetra- 
tion of  the  sphere.  We  shall  also  set 

The  flow  of  the  fluid  Is  due  to  the  effective  portion  of  the 
sphere  and  Is  approximately  due  to  a  lens  with  dimensions  B,, 
c.  In  place  of  B,c. 

The  dimenslonless  virtual  mass  m  =  — M —  is  given  by 

mj.(b^)  of  the  lens  flow*  Thus  the  dimenslonless  virtual  mass 

and  the  drag  coefficient  ^  are 

db 

(6.1) 

^^dm^(b^)     dbj^ 

db    6b^         '      db 

dmj.(b, ) 
where  m- (b, )  and  — ^- are  determined  from  the  lens  flow 

with  the  argument  b,» 

It  remains  then  to  determine  b,  as  a  function  of  b.  Set 

(6.2)  b^  =  b  .  w(b) 

(or  equlvalently,  B,  =  B  .w(b)  )  where  w(b)  Is  a  function  of 
b  =  B/A.  The  quantity  w(b)  will  be  called  the  "wetting  factor." 
In  this  section  we  shall  show  theoretically  that  Initially, 
for  b  =  0,  the  wetting  factor  w(0)  has  the  value  3/2.  The  exact 
dependence  of  w(b)  on  b  Is  extremely  difficult  to  determine 
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theoretically,  but  for  small  values  of  b  the  theory  shows 
that  w(b)  behaves  like  -  -  constant  ./b^   Consequently, 
the  variation  of  w(b)  with  b  was  determined  from  experimental 
me  a  sur  erne  n  t  s  • 

A  blackened  steel  sphere  of  diameter  2.86  cm  was  dropped 
into  a  tank  of  water  from  a  height  of  16.2  cm  and  flash  photo- 
graphs were  taken  with  the  sphere  at  various  stages  of  Immersion. 
The  photographs  In  Plates  2  and  3  are  a  selection  showing  the 
sequences  of  events.   Plate  4  shows  an  enlarged  view  of  the 
splash  and  Illustrates  how  the  measurements  were  made.   The 
straight  line  represents  the  Initial  surface  level,  the  posi- 
tion of  which  was  determined  by  a  pair  of  inclined  needles  and 
their  reflections,  as  indicated  in  the  photograph.   The  boundary 
of  the  sphere  has  been  con^leted  by  a  white  line.   The  exact 
line  of  demarcation  between  sheath  and  base  is  somewhat  arbitra- 
ry but  estimates  can  be  made  with  a  fair  degree  of  accuracy  and 
consistency.   The  upper  pair  of  short  horizontal  lines  shows  our 
choice  of  demarcation  between  sheath  and  base.   Below  this  line 
the  risen  water  is  thick  while  above  it  the  water  thins  out  very 
rapidly.  For  this  particular  photograph,  the  values  of  b  and 
w(b)  are: 

b  =  .27 

w(b)  =  1.30   . 

Similar  measurements  were  made  on  a  number   of   enlarged 
photographs   showing  various   depths   of  penetration.      The  quantity 
w(b)  was  then  plotted  against  /F  with  the  results   indicated  in 
Graph  8(a).      Also  plotted  in  Graph  8(a)   is   an  experimental  point 
obtained  by  Watanabe  while  measuring  impact  forces.      This    last 
point   is   discussed  in  Section  9. 


*  The   experiments  were  performed  at   the  AMG-NYTJ  Laboratory  by 
M.   Shamos. 
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The   experimental  points   shown  In  Graph  8(a)    fall  on 
a  reasonably  well-defined  straight  line  passing  through  the 
theoretical  value  w(0)    =  3/2.      The   equation  of  this    line   Is 


(6.3) 


w(b)    =  1.5     -     .4  /b" 


In  Graph  8(b),  w(b)  Is  drawn  as  a  function  of  b.  Actually, 
In  Graph  8(a),  the  experimental  points  are  all  slightly  above 
the  straight  line  (6«3)  near  /S  =1.  A  more  accurate  curve 
fitting  the  data  is  the  parabola 

w(b)  =  1.5  -  .4/b'  +  .065  b  . 


But  this  will  make  a  noticeable  difference  only  near  b  =  1,  and 
we  shall  rather  use  the  straight  line  (6.3). 

It  remains  to  derive  the  theoretical  value  w(0)  =  3/2  • 
We  shall  use  a  method  similar  to  one  applied  by  Wagner  [16] 
for  the  case  where  the  entering  body  is  a  nearly  flat  wedge 
(two  dimensional).  In  this  case,  the  wetting  factor  turns  out 


to  be  4i  =  1.57< 


A  corresponding  calculation  for  a  very  flat 
4 
cone  gives  a  wetting  factor  —  =  1.27. 

Tf 

Consider  the  wetting  correction  for  small  b,  when  the 
sphere  has  entered  only  a  small  distance  into  the  fluid.  The 
flow  of  the  fluid  is  approximately  due  to  the  spherical  segment 
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QQ«   Since  B  and  B^  are  supposed  small  the  flow  due  to  this 
spherical  segment  can  be  approximated  by  that  due  to  a  cir- 
cular disc  of  radius  c,»  These  approximations  are  accurate 
In  the  limiting  case  as  B  — c>  0,  as  can  be  seen  by  the  follow- 
ing argument.   Perform  a  similarity  transformation  with  magni- 
fication factor  —  .   The  limiting  situation  as  B  -»►  0  is  indi- 
cated by  the  diagram  below.   The  effective  portion  of  the  sphere 


in  contact  with  the  water  becomes  the  circular  disc  QQ,  while 
the  free  surface  becomes  the  plane  Q  oo  • 

The  theory  of  the  flow  due  to  a  circular  disc  is  classi- 
cal (  [8]^  pp.  130,  135,  or  [10],  pp.  456,  457).   The  potential 
function  1//  describing  the  flow  is  given  in  [8],  p.  135.   In 
particular  the  value  of  y^   on  the  lower  siirface  of  the  disc  is 


(6.4) 


r 


2B  r^ 
irv^i  - 


where  r  is  the  distance  from  the  center  and  c,  is  the  radius 
of  the  disc.   The  value  of  the  upward  velocity  -  -^    on  the 
free  surface  of  the  fluid,  level  with  the  disc,  is  (  [8], p.  130) 


(6.5) 


|2f  «2B[ 


lA 


/l-(v^)' 


.in  -^   ] 
r 
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The  motion  of  the  surface  of  the  fluid  can  now  be 
followed.  At  any  point  of  the  free  surface  at  distance 
r  from  the  axis  of  symmetry,  let  •)]  (r)  be  the  height  of  the 
free  surface  at  the  time  t  above  its  original  level  (the 
dependence  of  ')7(r)  on  the  time  t  is  understood.  Asstame  that 
the  upward  velocity  of  the  surface  at  the  Instant  t  is  that 
induced  by  the  disc  of  radius  c^j^.  Prom  (6.5),  we  have 

(^•6)  ^  (r)  =  BF  (1L) 

where 

(6.7 )  p(^)  =  I  [  ,^  -  arc  sin  I  ]  . 


^i-t 


The  quantity  c^  is  a  function  of  B, ,  the  relation  between 
them  depending  on  the  shape  of  the  Incoming  object.  For  a 
sphere  of  radius  A  we  have 

(6.8)  c^  =  2AB^   -  B^  • 

For  small  B,  (6.2)  can  be  written  as 
B^  =  B  .  w(0)   , 
and  (6«8)  then  is 

(6.9)  c^  =  2AB^  =  2A  w(0)  .  B   • 

^^^  .   c,c, 

B  =  -3^-i-  , 
Aw(0) 

and  (6.6)  becomes 

r*^i      c 

(6.10)  il(r)  =  — ^—     P(-l)cTdc,  . 

Aw(0)/      r    -L  J- 
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(The  upper  limit  c^  of  Integration  should  be  smaller  than 
r»   The  fluid  particle  initially  at  a  distance  r  rises  until 
It  hits  the  sphere,  after  which  it  enters  into  the  thin  sheath). 
In  this  integration  r  is  kept  fixed.   Set 

(6.11)  -1  =  ^ 

r 

where  ^  -  1,   Then  (6.10)  becomes 

07(r)  =-E fn^)idi     , 

Aw(0)  J^ 

or,  using  (6.11)  and  (6.9), 

(6.12)  -AM   =  B  .  2   /^p(n  i   d^. 

Equation  (6.12)  represents  the  equation  of  the  surface 
as  long  as  r  -  o-,.  In  particular,  at  the  position  r  =  c-,,  or 
^  =  1,  the  height  of  the  surface  is 

B^  -  B  =  B  .  [w(0)  -  1]  . 

Equation   (6.12)   yields 

(6.13)  w(0)    -   1  =  zjnl)   l    dU 

Substitution  of  (6.7)  into  (6.13)  gives  the  value  of 
the  wetting  factor  w(0)  at  b  =  0: 

(6.14)  w(0)  =  I  . 

Evaluation  of   (6.12)   gives  the  equation  of  the  surface  of  the 
water 

(6.15)  7f(r)   =  B   .  i  [    (y2   -  2)   arc  sin  ^    -  | /l   -  ^^3    . 
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Equation  (5. 14)  Is  plotted  below  for  the  case  when  "b  =  .01. 
(Of  course,  this  method  cannot  be  expected  to  produce  the 
thin  sheath* ) 


Section  7.   The  Free  Surface  Correetion 

The  free  surface  correction  arises  from  the  term 

(7.1)  /'//Vdxdy 

S 
occurring  in  (4«3)  and  the  fact  that  we  shall  no  longer  neglect 
squares  of  velocities  on  S»   The  exact  boundary  condition  is 
then 


(7.2) 


8£  _  1 


at 


(gradT/r)' 


in  place  of  y  =  0   (see   (5.1)   and  the  accompanying  discussion). 
Suppose  we  follow  a  particle  of  fluid  on  the  free   surface 
and  determine  the  rate  at  which  the  value  of   V'   changes  for 
that  particle.      This  requires  use  of  the  "particle  derivative" 
23^  where 


Dt 


Dt      at     ®      ^ 


grad  Y  . 
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In  view  of  (7»2),  we  have 

(7.3)  m^'k  ^^^^y^^^ 

on  the  free  siirface  S«  The  value  of  V^  on  S  thus  constantly 
decreases  and  Is  negative  at  any  Instant  after  the  Initial 
contact  of  the  sphere  with  the  fluid.   This  shows  that  the 
surface  correction  (7#1)  Is  negative. 

The  correction  to  the  Impact  force  P  arises  from  the 
time  derivative  of  (7.1).   Let  the  height  of  the  siirface 
above  the  Initial  level  be 

'>?  =  '>?  (r,t). 

Denote  the  surface  correction  to  the  Impact  force  Dy  P^  ,  so 
that  from  (4.3), 

(..4)   P,  =,//|ra.a.=,//{f.,|£}a..y  . 

5  S 

(There   Is   no   contribution  from  the   changing  domain  of   Integra- 
tion    since   this   is   cancelled  by  the   corresponding   term  from 

4^   (/°//yclx<iy)    •)      Substitution  of   (7.2)    gives 

w 

(7.5)  P3  «/'jJ|i(gradr)^  *  "^  |f  |  dxdy     . 

Now  asstime  as  a  reasonable  approximation  that  the 

velocity  on  the  surface  S  is  that  Induced  by  the  lens  flow, 

used  as  a  first  approximation  in  Section  5.   The  velocity  of 

the  surface  is  then  con^^letely  vertical,  so  that  V  =  -  rr^ 

oz 
and  (7.5)  becomes 


<^-^'  ^s  =  -i^//(||)'^^T  =  -|B 


^  )   dxdy 
oz 
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The  velocity 


32 


is  Infinite  at  the  corner  of  the  lens. 


but  the  integral  (7.6)  nevertheless  converges.  However, 
the  Integral  (7.6)  is  a  bit  cumbersome  to  evaluate  for  a 
lens  flow;  we  shall  Instead  approximate  the  flow  at  each 
Instant  by  that  Induced  by  an  oblate  spheroid  circumscrib- 
ing the  lens.   (In  the  previous  report  [5],  this  oblate 
spheroid  was  shown  to  be  a  reasonable  approximation  to  the 
lens. ) 

The  flow  about  an  ellopsold  is  well  known  (see  [10], 
pp.  456-7)  and  [8],  p.  132  ff.  ). 


The  velocity  nj    on  the  surface  is 


r/  a  B 


2-f(l) 


where 


'( t  )  =  -^  f-p==L===2  -   ,  ^   ..  arc  sin  (  ]/l^  |)\ 


.  C,         B- 

and  ^  =  —  ,  V  =  — -  ,  The  result  of  substituting  this  into 


(7.6)    is 
(7.7) 


Pa=-   4/oB       .       Tf    A*^    .    Kg 


S~"  2/° 
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where  r,  f         -,        ^  '-j^' 

^4|log(l)-l-« 


A^  )  c  arccoaj?  ^-i   2        2j 
/l  -  V' 

Kg  Is  drawn  as  a  function  of  b  =  B/A  In  Graph  ?•   This  correc- 
tion Is  not  to  be  considered  as  accxirate  for  very  small  b,  but 
should  be  a  reasonable  estimate  over  most  of  the  range  In  which 
we  are  interested. 

Section  8.  The  Second  Appro3tlmatloit 

The  results  of  Sections  5,6,7  will  now  be  collected  to- 
gether. Set 

W 

SO  that  M  Is  the  contribution  to  the  virtual  mass  (4.5)  due 
to  the  wetted  region  W  alone.   Then  eqiiation  (4.3)  can  be 
written  in  the  form 


P  =  «    (V)  *  f. 


where  P  Is  given  by  (7.4),  or 


s 

n    dM^  .2 
"dB 


V+^S  *P. 


Using  MB  =  -  P  and  (7.7),  we  have 


or 


(8.1)        P=iou2.wA^.  ^ 5 


2'  o 
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M  3 

In  this  formula,  the  denominator  can  be  approximated  by  (1  +  —  )  • 

Comparing  formula  (8.1),  the  denominator  being  replaced 
hy  (1  +  M  )3^  with  formulas  (3.6),  (3.7),  we  see  that 

Mo 


db   db     a 


where  m  is  the  final  estimate  of  the  virtual  mass,  and 

=V  =  =»L^^1^   • 

Thus  m  may  be  obtained  by  integration,  with  the  result 

ob 
(8.3)  m  =  mj^(b^)  -/   KgC 


.db 


This  is  drawn  in  Graph  3  and  conipared  with  the  first  approxima- 
tion iBT(b)  arising  from  the  lens  flow. 

The  impact  drag  coefficient  Is  calculated  from  formula 

(3.7)  using  the  final  estimate  of  m.  The  numerator  ^  can 

db 
be  written,  by  (8.2),  in  the  form 

^„   dmT(hT)   db, 

db     db-,     db     s 

dmT(b,  ) 
The  quantity  — = — =-  may  be  read  from  Oraph  2,  which  is  a  repeti- 
db^ 

^  db, 

of  the  result  obtained  in  [15];  ■-— =  is  obtaiiied  from  formulas 

db  Ajf, 

(6.2),  (6.3);  and  K.  is  read  from  Graph  7.  The  result  for  ^ 
a  •  aO 

is  drawn  in  Graph  1  (the  case  r  =  oo  ),  This  is  exactly  the 

impact  drag  coefficient  CU  in  case  the  mass  M^  of  the  Incoming 

projectile  is  large  compared  to  the  mass  of  the  water  displaced 

by  the  hemispherical  nose. 
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The  lirpact  drag  coefficient  Cp  In  the  general  case 

of  finite  6"  may  be  calculated  from  forimila  (3»9)  and  the 

graphs  of  m  and  — •   The  results  corresponding  to  r  =  cxs  ^ 

db 
1,  #5,  •25  are  shown  In  Graph  1. 

We  observe  that  the  maximum  value  of  C_  Is  about  1 
(unless  (T  Is  very  email),  and  that  this  maximum  occurs  for 
b  =  .IS  approximately. 

In  the  neighborhood  of  b  =  0  the  expansion  of  Cp  can 
be  obtained  from  (8.4),  using  (5»3)  or  (5»4).   The  result  is 

(8.5)  C  =  6.62  b^/^  -  agb  +  ... 

where  ag  is  a  coefficient  which  is  not  accurately  given  by 
(8»4)   and  Its   determination  is   difficult.      The   coefficient 


6.62   is  1  .    3^/2. 

IT 

In  Graph  2,  the  final  Cp  is  compared  with  the  Cp  ob- 
tained from  the  first  approximation  of  the  preceding  report [15]. 

Section  9.  Coigparison  with  Experiment 

Careful  experiments  on  the  impact  of  cones  and  spheres 
have  been  conducted  by  S.  Watanabe  (see  [17].  He  used  a  thin 
spherical  segment,  with  dimensions  as  indicated  in  the  diagram, 
on  top  of  which  were  placed  a  weight  and  piezoelectric  gauges 
connected  to  an  oscillograph.  The  system  was  dropped  from 
heights  of  20  cm  to  50  cm  from  the  water  surface,  and  the 
oscillogram  record  gave  directly  the  impact  force  as  a  function 
of  the  time.   The  Watanabe  results  are  exhibited  in  Graph  8(a), 

\ 
h  =  .84  cm.        \ 
Mass  M^  =1.645  kgra.       ^ 

r  =  — 2—  =  .116  \ 

. ^ V-^ 

\h  =  >'■/  cm . 
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takon  from  his  paper,  the  vertical  axla  R  being  the  Impact 
force  ineasuped  In  kilograms  of  force  and  H  being  the  height 
from  which  the  spherical  segment  was  dropped* 

These  graphs  are  converted  Into  graphs  of  C  as  follows: 


Jfi. 


|/oTj2    .wA^        i^,   2gH  .itA^         lyjA^H 


.    depth  of  Immersion 


where  the  radius  A  of  the  sphere  Is  15  cm.   The  results  are 
exhibited  In  Graph  8(b). 

Theoretically  one  should  expect  the  formula  (3»9)  for 
Cp  to  be  valid*  Expanding  m  and  ^  In  powers  of  b, 

m  =  I  a^bV^  -  ^b^  +  ... 

(9.1) 

'  dm  ^  „  .1/2       ^  y,     ^ 
--  =  a^b  /   -  agb  +  ... 


where  a,,  a^  are  constants  (from  (8.3)  the  theoretical  estimate 

3m  3        9m 
of  a.  Is  6.62.   The  denominator  (1  +  -=)  «  1  +  —  near  b  =  0, 

1  8(r        8(r 

80  that 

a^bV2  .  agb 
(9.2)  C 


4    <r      16   <r 


i  =  .^.V2[i  ||,..^Mi.^^] 
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Sotting 

(9.3) 

w«  have 
(9.4) 


,,.V3ti.||l, 


16  (T  ^ 


% 


Igx  «  1 


The  quantities  y  and  x  are  plotted  in  Graph  9(a) 
using  Watanabe^s  experimental  points.   They  show  an  excellent 
agreement  with  the  linear  relation  (9»4),  and  the  best  straight 
lines  are  likewise  drawn  in  Graph  9(a),  the  resulting  values 
of  a,  and  ag  being: 


H(cm) 

20 

30 

40 

50 

*1 

6.96 

6.92 

7.17 

7.13 

*2 

23.5 

25.9 

26.6 

27.3 

The  values  of  a,  show  good  agreement  with  the  expected  theoretl" 
cal  value  6.62. 

Selecting  in  equation  (9.2) 


(9.5) 


6.62 


and  dropping  the  denominator  completely,  we  have 

6.62  b-*-/^  -  c_ 


(9.6) 
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6.62  bV2  _  Q 

The  quantity  r ^  is  plotted  against  b  in  Graph  9(b), 

the  values  of  b  less  than  .01  being  excluded  since  any  experi- 
mental error  in  c  will  then  produce  a  large  error  in  ag. 

The  result  shows  that  the  quantity  (9»6)  is  accurately 
constant  and  its  best  value  from  the  experimental  data  is  23. 5. 
If  the  denominator  in  (9»2)  is  included,  a  correction  to  the 
value  of  ag  is  obtained  which  amounts  to  -  1»5,   Thus,  the  ex- 
perimental value  of  ag  is 

(9.7]  ag  =  22  . 

Thus  the  theoretical  result  (9.5)  is  verified  completely 
by  Watanabe»s  experiment.   In  particular,  the  value  ^  for  the 
wetting  correction  is  verified. 

Another  more  direct  check  of  the  wetting  correction  of 
9   is  provided  by  Watanabe  in  his  observation  of  the  depth  of 
immersion  at  which  the  resistance  curve  dropped  sharply.   The 
depth,  in  Table  IV  of  Watanabe's  paper  [17],  is  given  as  .SSScm 
on  the  average.  This  is  clearly  due  to  the  water  overflowing 
the  spherical  segment.  The  wetting  factor  is  therefore  approxi- 
mately 

w(b)  =  -1^  «  1.44 
.585 


which  is  slightly  less  than  3/2.  The  corresponding  value  of 
b  is 

.  _  .585  _  -.,Q 
b  =  — -—  =  .039  . 
15 

This  experimental  point  is  plotted  in  Graph  6(a). 

In  the  experiments  by  Blundell  [3],  accelerometers  with 
low  natural  frequencies  of  100-400  cycles/sec  and  low  damping 
were  used.   The  responses  of  the  instruments  were  oscillatory, 
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80  that  only  the  general  shape  of  the  Impact  ciirve  could 
he  estahllshed.   The  dimenalona  of  the  dropped  sphere  are 
as  indicated* 

(A)  M^  =  351b,  Uq  =  11,4  ft/aec,  (T  =  .23 

(B)  M^  =  321b,  U^j  =  11,4  ft/sec,  (T  =  ,21 

(C)  M^  =  321b,  U^  =  8  ft/sec   ,  0"  =  .21 

(D)  M^  =  501b,  U^  =  8  ft/sec   ,  (T  =  .33 

The  experimental  ciirvea,  obtained  from  Plgtires  1,  2  of  [3] are 
converted  into  C  ciirves  in  Graph  11(a),  and  the  corresponding 
theoretical  curves  are  drawn  in  Graph  11(b).  An  idea  of  the 
experimental  error  can  be  obtained  by  comparing  (B)  and  (C), 
which  should  be  the  same,  and  by  noting  that  (D)  should  give 
the  highest  Cp  since  or    is  the  largest.  The  agreement  with 
theory  is  satisfactory,  and  curve  (D)  even  agrees  numerically 
with  theory. 

Richardson  [13]  used  a  hemispherical  shell  with  the 
same  dimensions  as  Blundell,  but  with  a  mass  of  24  lb  (  <r  =  .159) 
dropped  from  a  height  of  2  ft.  He  used  an  instrument  which 
worked  on  the  changes  in  capacity  of  a  "breathing"  condenser. 
His  experimental  graph.  Figure  2b  of  [13],  is  converted  into 
the  dlmenslonleas  form  of  a  Cp  c\irve  and  con5)ared  with  the 
theoretical  curve  in  Graph  12.  The  agreement  is  satisfacrory, 
especially  on  the  position  where  the  maxiraim  Cp  is  attained. 

In  all  these  experiments,  the  value  of  (T  is  rather  small. 
It  is  desirable  to  perform  further  experiments  with  a  larger  <r  , 
using  accurate  instruments  (such  as  piezoelectric  gauges). 
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Appendix 

Section  10»  Comparison  of  Energy  and  Momentum 

In  Sections  3  and  4  the  force  of  Impact  was  derived 
from  momentum  considerations,  and  we  found  that  the  force 
on  the  entering  sphere  Is  equal  to  the  time  rate  of  change 
of  the  resultant  Impulsive  pressure  force  acting  on  the 
boundary  of  the  fluid  (formula  (4,3)  )•  We  may,  however, 
derive  the  force  of  Impact  from  energy  considerations. 

Let  i  M  "  be  the  kinetic  energy  of  the  fluid  at  time 
t,  where  M^  Is  the  "energy"  virtual  mass.  By  the  principle 
of  conservation  of  energy,  we  have 

(10.1)  I  (M3  +  M^)  TJ2  =  I  M^TJ^  ,  or  IJ^  .  _2_  . 

1  +  -^ 
«o 

Differentiating  with  respect  to  t  and  using  (3*2 ),  we  obtain 

dM  dM^ 

(10.2)  p=i-li-n2  =  i_3I j^ 

l+_±       (1  +  -^) 
Mo  Mo 

by  (10.1). 

Formulas  (10.1),  (10.2)  should  be  compared  with  the 
corresponding  formulas  (3.3),  (3.4)  based  on  momenttm  con- 
siderations. An  Immediate  comparison  of  (10.1)  and  (3«3) 
shows  that 

(10.3)  1  +-S=  (1  +5  )2  or  M.  =2M(l+i4) 

"o       "o  o 
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where  M  Is  the  "momentim"  virtual  mass  which  we  have  been 
using  throughout*   The  same  relation  Is  obtained  by  compar- 
ing (10.2)  and  (3«4).  For  incoming  projectiles  of  large 
mass  compared  to  the  mass  of  the  water  displaced  by  the 
hemispherical  nose,  ^  M/M^  Is  small  compared  to  1,  and  we 
have  approximately 

(10»4)  M^  =  2M  . 

Thus,  for  impact,  M^  is  different  from  M. 

In  order  to  obtain  an  expression  for  M  ,  we  observe 
that  the  total  kinetic  energy  of  the  fluid  is 

by  Green's  theorem,  where  the  differentiation  is  with  respect 

to  the  normal  pointing  into  the  fluid*   Thus  the  kinetic 

energy  is  equal  to  the  work  done  by  the  impulsive  forces  acting 

on  the  boundary  of  the  fluid.  Since  -  rr^  =   U  cos<x  on  W  (see 

diagram  on  page  17),  we  see  that  -  ~  dS  =  U  dxdy.  Substitu- 

on 
ting  from  (4.4 ),  we  therefore  obtain  for  M  the  formula 

(10.5)       Mg  ^/°ff^  ^^  *  /°ff^   ^'  M^'^   • 
w  s 

Comparing  the  right  hand  sides  of  (10.5  and  (4«6),we 
see  that  they  differ  only  in  their  second  terms,  the  free  surface 
terns.   If  the  same  first  approximation  is  made  as  previously, 
that  is  ^  =  0  on  the  surface,  we  would  obtain 

"e  =  rff^  ^^-y   =  M 

k^ 

*  It  is  easily  seen  (from  the  considerations  of  Section  4)  that 

the  integration  extended  over  a  large  spherical  membrane  tends 

to  zero  as  the  membrane  tends  to  infinity. 
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whereas  we  know  that  M^  =  2M  approximately.  The  complete 
e 

neglect  of  the  surface  terms  may  thus  cause  a  sizeable 
error  amounting  to  the  factor  2» 

We  shall  show  that  the  neglect  of  the  surface  term 
is  justifiable  when  momentum  Is  used  and  not  justifiable 
when  energy  Is  used.  The  reason  for  this  Is  that  the  rapid- 
ly rising  fluid  Immediately  surrounding  the  sphere  carries 
considerable  energy  but  not  much  momentum.  This  portion  of 
the  fluid  includes  the  thin  sheath  and  the  thick  base  of  the 
surface  rise  adjacent  to  the  sphere,  all  of  which  will  be 
grouped  together  and  called  the  splash.  We  may  thus  think 
of  the  splash  as  carrying  the  energy  which  would  normally  be 
dissipated  into  heat  in  inelastic  Impact. 

We  shall  show  that  the  percentage  of  energy  carried 
by  the  splash  is  much  larger  than  the  percentage  of  momentum 
carried  by  it,  at  least  for  the  depths  of  penetration  in 
which  we  are  interested.  This  is  plausible  on  intuitive 
grounds  since  the  mass  of  the  splash  is  concentrated  near  the 
sides  of  the  sphere  and  is  moving  outward  as  well  as  upward 
(see  Plate  5  ).  The  outward  or  horizontal  motion  of  the  splash 
contributes  to  the  kinetic  energy,  but  contributes  nothing  to 
the  momentum  because  of  the  rotational  symmetry  about  the  verti- 
eal  direction. 

At  each  instant,  the  additional  water  entering  the  splash 
is  shed  off  the  side  of  the  sphere  at  the  place  where  the  water 
level  intersects  the  sphere.  Let  u,  and  vu   be  the  horizontal 
and  vertical  con^jonents  of  the  splash  as  it  leaves  the  sphere, 
and  let  8  be  the  angle  Indicated  in  the  diagram.  The  velocity 
with  which  the  splash  is  moving  in  the  direction  normal  to  the 
sphere  is  TJ  cos  ©,  where  TJ  is  the  vertical  velocity  of  the 

*  This  remark  is  not  valid  in  the  case  of  oblique  entry. 

Because  of  the  asymmetry  there  is  then  a  sizable  horizontal 
component  of  the  momentum  of  the  splash* 
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sphere.     We  therefore  have   {see  the   diagram) 
Vi^  =  (u^  +  U)   cot  e 

and  30  the  magnitude  u  of  the  resultant  splash  velocity  Is 
given  by 

u  =  u^Jl  +  (1  +~;)^  cot^  ©     . 


Let  M  be  the  virtual  mass  of  the  water  defined  In 
Sections  3,  4.   By  (10.4),  M^  =  2M  and  so  the  total  kinetic 
energy  of  the  fluid  Is  ^  M^uS  =  MU^.  Hence,  If  mg  Is  the 
mass  of  water  entering  the  splash,  the  fraction  of  the  energy 
carried  by  the  splash  Is 

(10.6)   -£i|^=l^A)2[l*(l+2)2eot2e)  . 

The  total  momentum  of  the  splash  Is  directed  virtually  upward 
by  rotational  symmetry  and  is  equal  to  nig„%»   Dividing  by 
the  effective  momenttim  MTJ  of  the  fluid  we  obtain 


(10.7) 


MU     1 


CONPIEENTIAL 
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for  the  fraction  of  the  moment\jm  carried  by  the  splash. 
The  ratio  of   (10,6)    to   (10.7)   Is 

(10.8)  percentage  enerp  In  splash     ^ 

*  '  percentage  momentum  In  splash 


l^[i+  (i^i   )2eot29]   >—l— 
'^   ^  \  sec  9-1 


"v 


for  any  -^   •  The  following  table  exhibits  the  minimum  values 
of  the  ratio  (10.8)  for  various  angles  0: 


minimum  possible 
value  of  ratio  (10.8) 


oo  28.3   6.5   2.4    1     .35    0 


This  table  gives  the  minimum  value  of  the  ratio  (10.8)  as  the 

fluid  is  shed  off  into  the  splash.  The  ratio  (10.8)  for  the 

total  splash  is  some  sort  of  integral  of  ^ .   Thus,  in 

sec  9-1 
the  early  stages  of  the  entry  of  the  sphere  the  energy  loss 

In  the  splash  is  much  greater  than  the  momentum  loss.  For 

this  reason,  an  approximation  based  on  momentiun  is  more  acciu?ate 

than  one  based  on  energy. 

In  Sections  7,  8  of  the  present  report,  we  have  computed 

the  correction  to  the  virtual  mass  M  due  to  the  motion  of  the 

surface.  The  surface  correction  to  the  dimenslonless  virtual 

mass  m  is  Z'  K  cib  .  This  correction  exhibits  the  behavior 

Indicated  above.  At  the  beginning,  near  b  =  0,  the  correction 
is  relatively  small.  The  relative  magnitude  of  the  correction 
Increases  until  near  b  =  1  It  amounts  to  about  20  percent. 
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Prom  the  mathematical  side  we  can  see  why  the  surface 

correction  la  smaller  on  a  momentum  basis  than  on  an  energy 

basis.   The  surface  correction,  on  a  momentum  basis.  Is 

merely  />]] ^   dxdy.   On  an  energy  basis  the  surface  correction 

Is  by  (lo!5)  /offit   (-  |i  )dS  which  has  the  extra  factor  |i  . 

s  an  an 

The  surface  of  the  water  Immediately  adjacent  to  the  sphere 

Is  moving  up  rapidly,  so  that  -^  Is  large  there  a.nd  />ff^{'   •|^)dS 

/</. .  an  '  "v    on 

Is  larger  than  /<>JJ  (p  dxdy. 

Also,  since  the  surface  Is  rising,  -  — ^  Is  negative  (  n  Is 

the  Inward  pointing  normal)  and  so  /off^i-^^dS   Is  positive, 

remembering  that  ^  Is  negative  on  the  free  surface  .   Thus 

M  >  foff^   dxdy  <  M^   , 
w 

the  Inequalities  being  due  to  the  surface  terms. 
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Polished   serpt-ntirif  sphere  falling    14  cm.    into  water. 


This  plate  is  taken  from  A.  M.  Worthing ton,  "i 
Study  of  Splashes",  1908. 


Plate  1 
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SMOOTH  ENTRY 


Steel   sphere   of  diameter   2.86   cm.    dropped 
from   a    height   of   16.2   cm. 


Plate   2 
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SMOOTH  ENTRY    (Continued) 


Steel   sphere   of   diameter   2.86   cm.   dropped 
from  a    height  of   16.2   cm. 


Plate   3 
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These  pictures  were  taken  by  Dr.  Hurley 
of  the  Morris  Dan  Group  at  California  Institute  of  Technology. 


Plate  5 
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ROUGH  ENTITY 


Steel   ball,    5/8"   diameter-,    velocity   19.4  ft/sec 
Pressure    above    liquid   =  20  cm,    mercury. 


These    pictures   were    taken  by  R.   M.   Davies 
of   the  Engineering   Laboratory,    Cambridge,   England 


Plate    6 
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HIGH-SPESD  ENTRY   SHOV/ING   COMPfffiSSIBILITZ  EFFECT:: 


U     =  8.84  X   10      cm/s 
o  / 


0.159   cm.      M     =0.6  grains. 


This   photograph  was    taken   by  E.    Newtor   Harvey 
of   Princeton   University. 
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